Direct Method - Quadratic 


Example 2 


The upward velocity of a rocket is given as a function of time in 
Table 2. 


Table 2 Velocity as a function of time. 


t (s) v(t) (m/s) 
0 0 

10 227.04 

15 362.78 

20 517.35 
225 602.97 

30 901.67 


Determine the value of the velocity at t=16 seconds using the 
direct method of interpolation and a second order polynomial. 


Solution 


For second order polynomial interpolation (also called quadratic 
interpolation), the velocity is given by 
v(t) =a,+atrt a,t’ 
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Figure 4 Quadratic interpolation. 


Since we want to find the velocity at t=16, and we are using a 
second order polynomial, we need to choose the three data points 
that are closest to t=16 that also bracket t=16 to evaluate it. 
The three points are ¢, =10, £ =15, and ¢, = 20. 
Then 

t, =10, v(t,)= 227.04 

t, =15, v(t,)=362.78 

t, =20, v(t, )=517.35 


gives 
v(10) = a, +a, (10)+a, (107 = 227.04 
v(15) =a, +a, (15)+a, (15) = 362.78 
v(20) = a, +a, (20)+a, (207 = 517.35 


Writing the three equations in matrix form, we have 
1 10 100| a, 227.04 


1 15 225]) a, |=| 362.78 
1 20 400}a,| | 517.35 


Solving the above three equations gives 
a, = 12.05 
a, =17.733 
a, = 0.3766 


Hence 
v(t)=12.05 + 17.733¢ + 0.376607, 10<t<20 


At t=16, 
v(16) = 12.05 +17.733(16) + 0.3766(16)° 
= 392.19 m/s 


The absolute relative approximate error obtained between the 
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results from the first and second order polynomial is 
7 Pex — 393.70 
-| 392.19 

= 0.38410% 
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